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Abstract — An asymptotically optimal trellis-coded modulation 
(TCM) encoder requires the joint design of the encoder and 
the binary labeling of the constellation. Since analytical ap- 
proaches are unknown, the only available solution is to perform 
an exhaustive search over the encoder and the labeling. For 
large constellation sizes and/or large number of encoder states, 
however, an exhaustive search is unfeasible. Traditional TCM 
designs overcome this problem by using a labeling that follows 
the set-partitioning principle and by performing an exhaustive 
search over the encoders. In this paper we study binary labelings 
for TCM and show how they can be grouped into classes, which 
considerably reduces the search space in a joint design. For 8-ary 
constellations, the number of different binary labelings that must 
be tested is reduced from 8! = 40320 to 240. For the particular 
case of an 8-ary pulse amplitude modulation constellation, this 
number is further reduced to 120 and for 8-ary phase shift keying 
to only 30. An algorithm to generate one labeling in each class 
is also introduced. Asymptotically optimal TCM encoders are 
tabulated which are up to 0.3 dB better than the previously best 
known encoders. 

Index Terms — Binary reflected Gray code, bit-interleaved 
coded modulation, coded modulation, convolutional encoder, 
performance bounds, set-partitioning, trellis-coded modulation, 
Viterbi decoding. 

I. Introduction 

The first breakthrough in coding for the bandwidth-limited 
regime came with Ungerboeck's trellis-coded modulation 
(TCM) m-H in the early 80s where the concept of labeling 
by set-partitioning (SP) was introduced. TCM was quickly 
adopted in the modem standards in the early 90s and is 
a well studied topic 121, E Sec. 8.12], iH Ch. 18]. The 
next breakthrough in coded modulation (CM) came in 1992, 
when Zehavi introduced the so-called bit-interleaved coded 
modulation (BICM) (H, 191, usually referred to as a pragmatic 
approach for CM fW\- 

The design philosophies behind TCM and BICM for the 
additive white Gaussian noise (AWGN) channel are quite dif- 
ferent. Ungerboeck's scheme is constructed coupling together 
a convolutional encoder and a constellation labeled using the 
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SP principle. For constellations having certain symmetries, SP 
can be achieved by using the natural binary code (NBC) 121 
Fig. 4], f?, Fig. 3]. On the other hand, BICM is typically a 
concatenation of a convolutional encoder and a constellation 
labeled by the binary reflected Gray code (BRGC) Ol, 113 
through a bit-level interleaver The BRGC is often used in 
BICM because it maximizes the BICM generalized mutual 
information for medium and high signal-to-noise ratios l9] 
Sec. Ill], flT, Sec. IV]. In TCM, the selection of the con- 
volutional encoder is done so that the minimum Euclidean 
distance (MED) is maximized, while in BICM the encoders 
are the ones optimized for binary transmission. These encoders 
are typically the maximum free Hamming distance codes [7] 
Sec. 12.3] or the so-called optimum distance spectrum (ODS) 
first tabulated in |T4 Tables III-V] and 1151 Tables II-IV] and 
later extended in |16|. 

It was recently shown in ITtII that if the interleaver is 
removed in BICM, its performance over the AWGN channel 
is greatly improved. This was later analyzed in detail in ifTSl 
for a rate R = 1/2 encoder and a 4-ary pulse amplitude 
modulation (PAM) constellation, where the system in ITTl 
was called "BICM with trivial interleavers" (BICM-T) and 
recognized as a TCM transmitter used with a BICM receiver 
Moreover, BICM-T was shown to perform asymptotically as 
well as TCM (in terms of MED) [2, Table I] if properly 
chosen convolutional encoders are used fl8l Table III]. The 
transmitters in 12 Table I] and flS^ Table III] for the 8- 
state (memory v = 3) convolutional encode^^ are shown in 
Fig. [T] (a) and Fig. [l] (c), respectively. 

The authors in fl8l| failed to note that in fact the optimal 
TCM encoder found when analyzing BICM-T is equivalent 
means to to the one proposed by Ungerboeck 30 years ago. 
For a 4PAM constellation, Ungerboeck's SP mapper (NBC) 
can be generated using the BRGC mapper plus one binary 
addition (which we call transform) applied to its inputs, as 
shown in Fig. [lib)- If the transform is included in the mapper, 
the encoder in Fig. (Ha) is obtained, while if it is included 
in the convolutional encoder, the TCM encoder in Fig. (He) 
is obtained. This equivalence also applies to encoders with 
larger number of stated and simply reveals that for 4PAM, a 
TCM transceiver based on a BRGC mapper will have identical 
performance to Ungerboeck's TCM if the encoder is properly 
modified, where the modification is the application of a simple 
transform. The idea of applying a linear transformation to the 

'Throughout this paper, all polynomial generators are given in octal. 

-We use "equivalent" to denote two encoders with the same input-output 
relationship. This is formally defined in Sec. |lll] 

^This equivalence does not directly hold because 1 18 Table III] Hsts the 
encoders in lexicographic order and because for some values of v there are 
more than one encoder with identical performance. 
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Fig. 1. Three equivalent TCM encoders: (a) convolutional encoder with 
polynomial generators G = [13,4] and an SP mapper |2|; (c) convolutional 
encoder with polynomial generators G = [13, 17] and a BRGC mapper ( 18|. 
The encoder in (b) shows how a transformation based on a binary addition 
can be included in the mapper (to go from (b) to (a)) or in the encoder (to 
go from (b) to (c)). 



labeling/encoder was first introduced in |fT9l Fig. 6.5] (see also 
II2OJ and L2L Ch. 2]). 

TCM designs based on SP are considered heuristic ll22l 
Sec. 3.4], and thus, they do not necessarily lead to an optimal 
design 125, p. 68O]0 The problem of using non-SP labelings 
for TCM has been studied in EH Sec. 13.2.1], ^ Sec. 8.6], 
and 126 1 . TCM encoders using the BRGC were designed in 
II24I . by searching over convolutional encoders maximizing the 
MED. In IM Ch. 6] and EqI, a non-Gray non-SP labeling was 
used and TCM encoders with optimal spectrum were tabulated. 

In a related work, Wesel et al. introduced in [27 1 the concept 
of the edge profile (EP) of a labeling, and argued that in most 
cases, the EP can be used to find equivalent TCM encoders in 
terms of MED. The EP is also claimed to be a good indication 
of the quality of a labeling for TCM in ll27l Sec. I]; however, its 
optimality is not proven. Consequently, an exhaustive search 
over labelings with optimal EP does not necessarily lead to an 
optimal design. 

In summary, as clearly explained in |[27l Sec. I], traditional 
TCM designs either optimize the encoder for a constellation 
labeled using the SP principle, or simply connect a convolu- 
tional encoder designed for binary transmission with an ad- 
hoc binary labeling. It has been known for many years that 
optimal TCM encoders are obtained only by jointly designing 
the convolutional encoder and the labeling of a TCM encoder 
Q p. 966]. However, to the best of our knowledge, there are 

**Indeed, the results in [24 Tables 2-3], (19 Ch. 6] and (20| show the 
suboptimality of the SP principle in terms of the multiplicities associated 
with the events at MED. 



no works formally addressing this problem, and thus, optimal 
TCM encoders are yet to be found. 

In this paper, we address the joint design of the feedforward 
convolutional encoder and the labeUng for TCM. To this end, 
we show that binary labelings can be grouped into different 
classes that lead to equivalent TCM encoders. The classes are 
closely related to the Hadamard classes introduced in ||28| in 
the context of vector quantization. This classification allows 
us to formally prove that in any TCM encoder, the NBC can 
be replaced by many other labelings (including the BRGC) 
without causing any performance degradation, provided that 
the encoder is properly selected. This explains the asymptotic 
equivalence between BICM-T and TCM observed in ifTSl . 
Moreover, since the classification reduces the number of 
labelings that must be tested in an exhaustive search, we 
use it to tabulate optimal TCM encoders for 4-ary and 8-ary 
constellations. 

II. Preliminaries 

A. Notation Convention 

Throughout this paper, scalars are denoted by italic letters x, 
row vectors by boldface letters x = [xi , . . . , x^], temporal se- 
quences by underlined boldface letters x— [a;[l], . . . , a;[A's]], 
and matrices by capital boldface letters X where Xi,j repre- 
sents the entry of X at row i, column j. The transpose of 
a matrix/vector is denoted by [-J^. Matrices are sometimes 
expressed in the compact form X = [xi; X2; . . . ; xm], where 
Xi = [xi^i, . . . , Xi^N] is the ith row. Sets are denoted using cal- 
ligraphic letters C and the binary set is defined as B = {0, 1}. 
Binary addition is denoted by a © 6. 

Random variables are denoted by capital letters Y, the 
probability mass function (PMF) of the random variable Y 
by Pyiv), and the probability density function (PDF) of the 
random variable Y by pviu)- Conditional PDFs are denoted as 
PY\x{y\x)- The tail probability of a standard Gaussian random 
variable is denoted by Q{x) = 1"^ d^. 

B. TCM Encoder 

We consider the TCM encoder shown in Fig. |2] where a 
feedforward convolutional encoder of rate R ~ k/m is serially 
connected to a mapper $£, and the index L emphasizes the 
dependency of the mapper on the labeling (defined later). At 
each discrete time instant n, the information bits ii „, . . . , „ 
are fed to the convolutional encoder, which is fully determined 
by k shift registers and the way the input sequences are con- 
nected (through the registers) to its outputs. Closely following 
the notation of fT. Sec. 11.1], we denote the length of the 
pth shift register by i^p, with p = 1, . . . , fc, the memory of 
the convolutional encoder by = '^p^ number 

of states by 2". The connection between the input and output 
bits is defined by the binary representation of the convolutional 
encoder matrix 11291 eq. (11.6)] 
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Fig. 2. Generic TCM encoder under consideration: A feedforward convolu- 
tional encoder of rate R = k/m with 2" states serially concatenated with a 
memoryless m-bit mapper 
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where gf = [9p\, ■ ■ ■ , 9pl^+iV £ 
vector representing the connection between the pth input 
sequence and the Zth output sequence with I — 1, . . . , m. The 
coefficients g^^\^ . . . , ffp j^^+i are associated with the input bits 
ip^m ■ ■ ■ ,ip,n-vp, respectively, and G G i^{'^+k)xm^ Through- 
out this paper, we will show the vectors gp ^ defining G either 
in binary or octal notation. When shown in octal notation, 
represents the most significant bit (see Fig. [T]i. 

The convolutional encoder matrix ([T]i allows us to express 
the output of the convolutional encoder at time n, which we 
define as m„ = . . . , Um,n], as a function of (z/ + k) 

information bits, i.e., 

Un=jnG, (2) 

where j„ = . . . , i^^^] with ijf^ = [ip,„, . . . , ip^n-i^J are 
the information bits, and the matrix multiplication is in GF(2). 

The coded bits m„ are mapped to real A^-dimensional 
constellation symbols using the mapper : — > X, 
where X C is the constellation used for transmission, with 
A"! = M = 2™. We use a;[n] e A" to denote the transmitted 
symbols at time n and the matrix X — [xi;x2;-.-;xm] 
with Xq e and q = 1,...,M to denote the ordered 
constellation points. We assume that the symbols are equally 
likely and that the constellation X is normalized to unit energy, 
i.e., ^ Ex[||Xf ] = 1/Af W^f = 1- The resulting 

spectral efficiency of the system is k [bit/symbol]. 

The binary labeling of the qth symbol in X is denoted by 
Cq = [cq^i, . . . ,Cq^m] £ 'B™, whcrc Cqj Is the bit associated 
with the kh input of the mapper in Fig. |2] The labeling matrix 
is defined as i = [ci; C2; . . . ; cj\/], where Cq in L corresponds 
to the binary label of the symbol Xq in X. Throughout this 
paper, we will show the vectors Cq in L in either binary or 
integer notation. 



j^i^p+i is a column 



C. Binary Labelings for TCM 

The NBC of order m is defined as N„i = [ni; 922; ... ; tim] 



where n„ 



] G S™ is the base-2 representa- 
tion of the integer q—1 and nq.,„ is the least significant bit. The 
BRGC of order m is defined as S,„ = [61; 62; • ■ • ; bi\i] where 
bq = [bq^i,...,bq^rn] £ B"\ The bits of the BRGC can be 
generated from the NBC as bq^i = rig^i and bqj = nqj^i^BUqj 
for / — 2, . . . , m. Alternatively, we have 

. . . ® bqj.l © bq^l for 1 = 1,. 



5g4 



Example 1: The NBC and BRGC of order m = 3 are 
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where the pivots of the labeling matrices (defined in Sec. IIII-BI l 
are highlighted. 

To formally define the SP principle for a given constella- 
tion X and labeling L, we define A;([m,„+i_;, . . . , m,„]) = 

{Xq £ X : [Cg^m+i_i, . . . , Cg^m] = , . . . , Um] , q = 

1, . . . , A/} C X for Z = 1, . . . , TO — 1. Additionally, we define 
the minimum intra-ED at level I as 



Si = min 



X, 



Xi 



, TO 



1. 



(5) 



and the MED of the constellation as ^o- 



Definition 1 (Set-partitioning [2]): For a given constella- 
tion X, the labeling L is said to follow the SP principle if 

So < 61 < . . . < 6„i-i- 

Example 2: Consider an 8PSK constellation (formally de- 
fined in Sec. |V]i. It can be easily verified that if this constella- 
tion is labeled by the NBC in (|4]i, an SP-labeled constellation 
is obtained. Although the NBC is the most intuitive form for 
generating an SP labeling for 7\fPSK constellations, it is not 
unique. As an example, consider the the semi set-partitioning 
(SSP) labeling proposed in ||30l Fig. 2(c)] and the so-called 
modified set-partitioning (MSP) labeling ISTl Fig. 2(b)]: 
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, TO, or, m matrix notation. 



It can be shown that both labelings follow the SP principle in 
Definition [1] 

Example |2] shows that there are multiple labelings that 
follow the SP principle. It can be shown that this is also the 
case for A/PAM constellations, and that in this case, the NBC 
is also an SP labeling. 
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D. System Optimization and Search Problems 

For a given constellation X and memory v, a TCM encoder 
is fully defined by the convolutional encoder matrix G and the 
labeling of the constellation L, and thus, a TCM encoder is 
defined by the pair O = [G,L]. 

For given integers fc, m, and v, we define the convolutional 
encoder universe as the set Qk.m,i' of {v + k) x m binary 
matrices G|f| which result in a noncatastrophic feedforward 
encoder and equally likely symbols]^ We are also interested 
in the labeling universe, defined for a given integer m as the 
set Cm of all M X m binary matrices whose M rows are all 
distinct. 

To the best of our knowledge, there are no works addressing 
the problem of designing a TCM encoder by exhaustively 
searching over the labeling universe and the convolutional 
encoder universe. We believe the reason for this is that an 
exhaustive search over encoders and labelings is unfeasible 
||32l Sec. I]. For example, for 8-ary constellations, there are in 
general 8! = 40320 different binary labelings. In this paper, 
we show how a joint optimization over all G G Qk,m.iy and 
L G Cm can be restricted, without loss of generality, to a joint 
optimization over all G e Qk,m,v and a subset of Cm- 
Ill. Equivalent Labelings for TCM Encoders 

In this section, we show that binary labelings can be 
grouped into classes, and that all the labelings belonging to 
the same class lead to equivalent TCM encoders. This analysis 
is inspired by the one in ll28l . where the so-called Hadamard 
classes were used to solve a related search problem in source 
coding. 

A. Equivalent TCM Encoders 

The transmitted symbol at time n of a given TCM encoder 
Q — [G, L] can be expressed using (|2]i as 

a;[n] =$i(w„) = $i(j„G). (7) 

Definition 2: Two TCM encoders Q = [G, L] and 9 = 
[G, L] are said to be equivalent if they give the same output 
symbol for the same information bit sequence, i.e., if they 
fulfill $z,(jG) = $i(jG) for any j G 8"+^. 

Remark 1: Two equivalent TCM encoders have the same 
bit error rate (BER) and frame error rate (FER). 

From now on we use T^j to denote the set of all binary 
invertible m x m matrices. 

Lemma 1: $i,(c) = $j^(cT) where L = LT, for any two 
mappers and $2, that use the same constellation X, any 
T eTm, and any c £ B"\ 

Proof: Let Uq = [0, . . . , 0, 1, 0, . . . , 0] be a vector of 
length Af, where the one is in position q. Thus, Cq = UqL 
for q = 1,...,M. The mapping $£, satisfies by definition 

^Note that whenever G is given in its binary form, vi, . . . are also 
needed to interpret G correctly according to (T). 

^For some matrices G, the symbols x[n\ can be nonequally likely. This 
would induce nonequally likely symbols (signal shaping) which we do not 
consider in this work. 



^L{cq) = Xq for q — 1, . . . , M, or, making the dependency 
on L explicit, 

$i(c)=a;,, if UqL (8) 
for any c e B™. Similarly, for any c e B™, 

<^l{cT)^Xq, if CT = UqL 

~ Xq, if C = UqL, (9) 

where the last step follows because L ~ LT^^. Since the 
right-hand sides of ^ and ^ are equal, ^i{cT) — $i,(c) 
for all c e B". □ 

The following theorem is the main result of this paper. 

Theorem 1: For any G G Gk.m.ui L G Cm, and T G %n, 
the two TCM encoders 6 = [G,L] and 6 = [G,L] are 
equivalent, where L = LT and G = GT. 

Proof: For any j G B''+'', 'i>i{jG) = ^i{jGT) = 
^L{jG), where the last equality follows by Lemma [T] The 
theorem now follows using Definition |2] □ 

Theorem[T]shows that a full search over Qk.m.v and Cm will 
include many pairs of equivalent TCM encoders. Therefore, an 
optimal TCM encoder with given parameters can be found by 
searching over a subset of Qk,m.w and the whole set Cm or vice 
versa. In this paper, we choose the latter approach, searching 
over a subset of Cm- 

B. Matrix Factorization 

We briefly summarize here some matrix algebra. The fol- 
lowing definition of a reduced column echelon matrix comes 
from l33] pp. 183-184], adapted to the fact fliat we only 
consider binary labeling matrices L whose columns are all 
nonzero. The first nonzero element of the fcth column is called 
the fcth pivot of L. The pivots for TVs and B3 are highlighted 
in ©. 

Definition 3: A matrix L G B^'^^™ is a reduced column 
echelon matrix if the following two conditions are fulfilled: 

1) Every row with a pivot has all its other entries zero. 

2) The pivot in column / is located in a row below the pivot 
in column I + 1, for / = 1, . . . , rn — 1. 

The matrix N3 in Example [T] (or more generally 7V„i) is 
an example of a reduced column echelon matrix. On the other 
hand, B,„ is not a reduced column echelon matrix because it 
does not fulfill the first condition in Definition [3] 

The following theorem will be used to develop an efficient 
search algorithm in the next section. We refer the reader to 
II33I p. 187, Corollary 1] for a proof. From now on we use 
TZm to denote the set of all reduced column echelon binary 
matrices. 

Theorem 2: Any binary labeling L G Cm can be uniquely 
factorized as 

L = LrT, (10) 

where T G Tm and Ljx G TZm- 

Theorem |2] shows that all binary labeling matrices L can 
be uniquely generated by finding all the invertible matrices T 
(the set 7m) and all reduced column echelon matrices (the 
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TABLE I 

Number of classes (Mr = |7?.m|), their cardinality 

(Mt = \Tm\), AND THE TOTAL NUMBER OF LABELINGS (2'"!) FOR 
DIFFERENT VALUES OF m. 



m 


1 


2 


3 


4 


5 




6 




Mr 


2 


4 


240 


1.038 • 10^ 


2.632 • 


^q28 


6.294 ■ 


10^« 


Mt 


1 


6 


168 


20160 


9.999 


10® 


2.016 • 


10^" 


2""! 


2 


24 


40320 


2.092 ■ 10^^ 


2.631 ■ 




1.269 ■ 


10«« 



set n,n). In particular, we have ||34l eq. (1)], ||28] eq. (18)] 

Mt^ |7;„| =n(2"-2'~'), (11) 



1=1 



(12) 



In Table H] the values for Mr and Mt for 1 < m < 6 
are shown. In this table we also show the number of binary 
labeUngs (\Cm\ = 2"M = MrA/t), i.e., the number of 
matrices L in the labeling universe. 

A modified Hadamard class is defined as the set of matrices 
L that can be generated via ( fTOl i using the same reduced col- 
umn echelon matrix Xr. Note that these modified Hadamard 
classes are narrower than the regular Hadamard classes defined 
in li28l . each including M reduced column echelon matrices. 
There are thus A/r modified Hadamard classes, each with 
cardinality Al^. 

As a consequence of Theorems [T] and |2] the two TCM 
encoders [G^L] and [GT^^jXr] are equivalent for any 
G G Gk,m.i> and L e £,„, where £r and T are given by 
the factorization ( fTOl i. In other words, all nonequivalent TCM 
encoders can be generated using one member of each modified 
Hadamard class only, and thus, a joint optimization over all 
G G Gk,m,L' and L G can be reduced to an optimization 
over all G G Gk,m,v and L G TZm with no loss in performance. 
This means that the search space is reduced by at least a 
factor of Mt = Ml/M^. For example, for 8-ary constellations 
(m = 3), the total number of different binary labelings that 
must be tested is reduced from 8! = 40320 to 240. Moreover, 
as we will see in Sec. |V] this can be reduced even further if 
the constellation X possesses certain symmetries. 

C. Modified Full Linear Search Algorithm 

The problem of finding the set TZm of reduced column 
echelon matrices for a given m can be solved by using a 
modified version of the full linear search algorithm (FLSA) 
introduced in ll28l Sec. VIII]. We call this algorithm the mod- 
ified FLSA (MFLSA). The MFLSA generates one member of 
each modified Hadamard class, the one that corresponds to a 
reduced column echelon matrix Xr . Its pseudocode is shown 
in Algorithm 1. In this algorithm, the vector r — [ri, . . . , tm] 
denotes the integer representation of the rows of the matrix Lr 
where = Cg,„i+2cg,,„_i + . . . + 2''""icg,i for g = 1, . . . , M. 
The first labeling generated (line 1) is always the NBC. 
Then the algorithm proceeds by generating all permutations 
thereof, under the condition that no power of two (1,2,4,...) 



Algorithm 1 Modified full linear search algorithm (MFLSA) 
Input: The order m 

Output: Print the Mr different reduced column echelon vectors r 
1 

2 
3 
4 
5 
6 
7 



[rM , index • ' • i '^M- 



9: 

10: 
11 
12 
13 
14 
15 
16 
17 



r ^ [0,1,...,M-1] 
loop 
print r 

index -S— 

while rM — index do 

[r'index+l, ■ ■ ■ , T'm] 

index <— index + 1 
while index is a power of 2 do 

index •<— index + 1 
end while 

if index — M — 1 then 

Quit 
end if 
end while 

Find pointer such that rpointer = index 
Swap rpointer and 

r pointer+l 

end loop 



is preceded by a larger value. By Definition [3] this simple 
condition assures that only reduced column echelon matrices 
are generated. 

Example 3: For to = 2, the MFLSA returns the following 
reduced column echelon matrices: 



7^2 



00 




01 




01 




01 


01 




00 




10 




10 


10 




10 




00 




11 


1 1 




1 1 




1 1 




00 



(13) 



where the first element in TZ2 is the NBC defined in Sec. III-CI 
and again we highlighted the pivots of the matrices. The 6 
binary invertible matrices for m — 2 are 



T2 



01 




01 




10 




10 




1 1 




1 1 


10 




1 1 


1 


01 




1 1 




01 




10 



(14) 



Using Theorem |2] all the 24 binary labelings in £2 (see 
Table can be generated by multiplying the matrices in TZ2 
and 72- 

Example 4: For to = 3, the MFLSA generates the reduced 
column echelon matrices (given in integer notation) shown 
in Table |ll] (read row by row). The first column in the table 
corresponds to the output of the FLSA of ||28l . Columns two to 
eight show the additional matrices generated by the MFLSA, 
which differ from the first column only in the position of the 
all-zero codeword. In this table we also highlight the labelings 
generated by the MFLSA that at the same time have optimal 
EP |27| for SPAM and 8PSK (see Sec. [V]). 

Example 5: If we study the labelings in Example |2l we 
find that the SSP belongs to the first modified Hadamard class 
(Xr — N3) while the MSP belongs to a different class, i.e.. 



isSP — -^3 



100 
010 
101 



^MSP 



1 1 1 

010 
001 



(15) 



where lJ^ — [0,1,2,4,7,6,5,3] (in integer notation) is the 
233th labeling generated by the MFLSA (see Table This 
shows that the NBC does not span all the labelings that follow 
the SP principle. 
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TABLE II 

Output of the MFLSA for m = 3, read row by row. The labelings shown in boldface have optimal EP for SPAM (first four 

COLUMNS) AND FOR 8PSK (FIRST COLUMN). 



01234567 


10234567 


12034567 


12304567 


12340567 


12345067 


12345607 


12345670 


01243567 


10243567 


12043567 


12403567 


12430567 


12435067 


12435607 


12435670 


01245367 


10245367 


12045367 


12405367 


12450367 


12453067 


12453607 


1 2453670 


01245637 


10245637 


12045637 


12405637 


I 2450637 


12456037 


12456307 


1 2456370 


1 245673 


10245673 


I 2045673 


I 2405673 


I 2450673 


12456073 


12456703 


1 2456730 


01234657 


10234657 


12034657 


12304657 


12340657 


12346057 


12346507 


12346570 


01243657 


10243657 


12043657 


12403657 


12430657 


12436057 


12436507 


12436570 


01246357 


10246357 


12046357 


12406357 


12460357 


12463057 


12463507 


12463570 


01246537 


10246537 


12046537 


12406537 


12460537 


12465037 


12465307 


12465370 


01246573 


10246573 


12046573 


12406573 


12460573 


12465073 


12465703 


1 2465730 


01234675 


10234675 


I 2034675 


12304675 


12340675 


12346075 


12346705 


12346750 


01243675 


10243675 


12043675 


12403675 


12430675 


12436075 


12436705 


1 2436750 


01246375 


10246375 


12046375 


12406375 


12460375 


12463075 


12463705 


12463750 


01246735 


10246735 


12046735 


12406735 


12460735 


12467035 


12467305 


12467350 


01246753 


10246753 


12046753 


12406753 


12460753 


12467053 


12467503 


12467530 


01234576 


10234576 


12034576 


12304576 


12340576 


12345076 


12345706 


12345760 


01243576 


10243576 


12043576 


12403576 


12430576 


12435076 


12435706 


12435760 


01245376 


10245376 


12045376 


12405376 


12450376 


12453076 


12453706 


12453760 


01245736 


10245736 


12045736 


12405736 


12450736 


12457036 


12457306 


12457360 


01245763 


10245763 


12045763 


12405763 


12450763 


12457063 


12457603 


12457630 


01234756 


10234756 


12034756 


12304756 


12340756 


12347056 


12347506 


12347560 


01243756 


10243756 


12043756 


12403756 


12430756 


12437056 


12437506 


12437560 


01247356 


10247356 


12047356 


12407356 


12470356 


12473056 


12473506 


12473560 


01247536 


10247536 


12047536 


12407536 


12470536 


12475036 


12475306 


12475360 


01247563 


10247563 


12047563 


12407563 


12470563 


12475063 


12475603 


12475630 


01234765 


10234765 


12034765 


12304765 


12340765 


12347065 


12347605 


12347650 


01243765 


10243765 


12043765 


12403765 


12430765 


12437065 


12437605 


1 2437650 


01247365 


10247365 


12047365 


12407365 


12470365 


12473065 


12473605 


12473650 


01247635 


10247635 


12047635 


12407635 


12470635 


12476035 


12476305 


12476350 


01247653 


10247653 


12047653 


12407653 


12470653 


12476053 


12476503 


12476530 



D. NBC and BRGC 

Another way of interpreting the result in Theorem [T] is that 
for any TCM encoder = [G, £], a new equivalent TCM 
encoder can be generated using an encoder G = GT^ and 
a labeling L = LT^^ that belongs to the same modified 
Hadamard class as the original labeling L. One direct conse- 
quence of this result is that any TCM encoder using the NBC 
labeling N,n and a convolutional encoder G is equivalent to 
a TCM encoder using the BRGC Bm and a convolutional 
encoder GT with T given by (|3]l. This is formalized in the 
following theorem. 

Theorem 3: The BRGC and the NBC of any order m 
belong to the same modified Hadamard class. 

Proof: The BRGC and NBC ai-e related via B,„ = 
N„iT, with T given by (|3]l. The theorem now follows from 
Theorem |2] and the definition of a modified Hadamard class. 

□ 

Example 6: For the two TCM encoders in Fig. [T] the NBC 
and BRGC labelings are related via B2 = N2T, i.e., 



"00" 




"00" 


01 




01 


11 




10 


10 




11 



1 1 

01 



(16) 



Thus, the BRGC and the NBC of order to = 2 belong to the 
same modified Hadamard class, and convolutional encoders 
can be chosen to make the two resulting TCM encoders equiv- 
alent. This was illustrated in Fig.[Tl where the transform block 
corresponds to the transform matrix T = [1, 1;0, 1] = T~^. 
Since N2 = B2T^, the TCM encoders [Gjig^iy] , B2] and 



[G[i3 4] , N2] are equivalent, where 

T 



G 



[13,4] 



1011 
0100 



= G 



[13,17]- 



1 _ 


'1 11 


T 


'1 1" 




1111 




01 



Example |6] and Theorem [3] explain, in part, the results 
obtained in ifTsl . where it is shown that the encoders in ifTsl 
Table III] used with the BRGC perform asymptotically as well 
as Ungerboeck's TCmQ 

IV. Error Probability Analysis 

The results in Sec. |lll]are valid for any memoryless channel 
model and any receiver; however, from now on we focus on the 
AWGN channel and a maximum likelihood (ML) decoder In 
this section, we briefly review bounds on the error probability 
of TCM encoders under these constraints. These bounds will 
be used in Sec. lIV-Bl to define optimal TCM encoders. 

Since TCM encoders are in general not linear, the proba- 
bility of error depends on the transmitted sequence, i.e., it is 
not possible to make the assumption that the all-zero sequence 
was transmitted ||5] p. 101]. This constraint can be lifted if the 
TCM encoder is "regular" ||35] Lemma 2], "superlinear" ll32l 
Sec. II-D], "scrambled" flSl, or "uniform" |[33, IH Ch. 18]. 
However, regularity, superlinearity and uniformity do not hold 
for all constellation and labelingo and thus, we cannot use it 
in this paper. 

We consider a baseband-equivalent discrete-time real-valued 
multi-dimensional AWGN channel. The transmitted sequence 
of equally likely symbols is denoted by ^ = [a;[l], . . . , cc[A^s]] 

^The "in part" comes from the fact that the system studied in ( 18| uses a 
(suboptimal) BICM receiver 

^For 8PSK for example, there is in fact no binary labehng that gives a 
regular' TCM encoder [22 Sec. 3.3], 
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where x[n] G X is the iV-dimensional symbol transmitted 
at discrete time n and Ng is the block length. The received 
sequence of symbols is y = [2/[l], ■ • • tUWs]], where y[n] ~ 
x[n] + z[n] G K,^ is the received vector at time instant n. The 
channel noise z[n] G is an A^-dimensional vector with 
samples of independent and identically distributed (i.i.d.) ran- 
dom variables with zero mean and variance No/2 per dimen- 
sion. The signal-to-noise ratio (SNR) is defined as E^/No — 
1/Nq. The conditional transition PDF of the channel is given 
by PY\xiy\xq) = (Not:)'^ ex.p{-No^^\\y - a;,|P). 

A. Error Bounds 

Let be the set of all length-i! symbol sequences that start 
at an arbitrary time instant and encoder state. Let Xe{x) be 
the set of length-i' sequences x^x that start and end at the 
same encoder state as x G and where all the other £ ~ 1 
intermediate states are different. An error event occurs when 
the decoder chooses a sequence x G Xi{x) which is different 
from the transmitted sequence x. Using the union bound, the 
probability of an error event of an ML TCM decoder at a 
given time instant can be upper-bounded as |22, eq. (6.6)0 



(17) 



= 1 xeXe 



ix) 

where PEP(x, i) is the pairwise eiTor probability (PEP) 
and Pxix) is the probability that the encoder generates the 
sequence x. 

Assuming i.i.d. information bits, the probability of the 
sequence starting at a given state is 1/2". There are 2^ equally 
likely branches leaving each state of the trellis at each time 
instant, and thus, 

1 1 

2^2* 

The PEP depends only on the accumulated squared ED (SED) 
between x and x and can be shown to be 



Px(x) = r^. 



(18) 



PEP(^,i) Q 



2No 



El 

n=l 



(19) 



Let Aj^2 j denote the number of pairs x G Xg and x G 
Xe{x) at accumulated SED d^ = X]l=i II^I"-] ^ ^i^d 
let A^ ,i2 f denote the number of pairs at accumulated SED d^ 
generated by input sequences at Hamming distance w. Using 
(fT8]l-(fT9]l and the definition of Ad2 j, ([TtI i can be expressed 



as 



where 



d^ev 



'2m 



(20) 



^"i" - E ^^^'^".^ " E E ^^^^ 



2" 2*= 

= 1 w=l 



'All the bounds in this section are dependent on the TCM encoder ©. 
However, to alleviate the notation, we omit writing out © as an explicit 
argument. 



is the distance multiplicity of the TCM encoder. In (|20] i V 
is the set of all possible accumulated SEDs between any two 
sequences, i.e., all the values of d^ for which A^2 ^ 0. 

To obtain a bound on the BER, each error event must be 
weighted by the number of bits in error (w out of fc), i.e.. 



BER < ^B^2Q 



Id^E, 

2Nn 



where 



1 1 \ - W 



2" 2^'^ ^ k ' 

= 1 w = l 



w,d^,e 



(22) 



(23) 



is the bit multiplicity of the TCM encoder. 

Finally, to obtain a bound on the FER we generalize the 
bound presented in ||37| for convolutional codes to obtain 



FER <NsY^ Ad2Q 
d^ev 



Id^Es 
'2K 



(24) 



B. Optimum Distance Spectrum TCM Encoders 

In this section we define TCM encoders that are optimal 
for asymptotically high SNR. These definitions will be used 
in Sec. |V] to tabulate optimized TCM encoders for different 
configurations. 

We call the infinite set of triplets {d^, A^2 , Bj2 } the distance 
spectrum (DS) of a given TCM encoder O = [G, L], where 
d^ G T). We also define the ith SED of a given TCM encoder 
by df with i = 1,2,3,..., where d^j^-^ > df and is the 
minimum SED of the TCM encoder. These SEDs correspond 
to the ordered set of SEDs in V. Based on ( l22l i and ( l24l i we 
define an optimum DS-TCM (ODS-TCM) as follows. 

Definition 4: A TCM encoder 9 = [G,L] with DS 
{d^,A42,Bii2} is said to have a superior DS to another TCM 
encoder 6 = [G,L] with DS {d^, Aj2, -8^2} if one of the 
following conditions is fulfilled: 

1) dl > dl or 

2) dl = dl, Aj^2 < Aj2 and B^2^ < B^2, or 

3) there exist an integer I > 1 such that d| = df, A^2 — A^2 
and B^2 = Sj2 for i = 1, 2, . . . , Z — 1 and df > df or 
df = df, A^2 < ij2 and ^^2 < B^2. 

Definition 5: For a given constellation X and memory 
the TCM encoder 9 = [G, L] is said to be an ODS-TCM 
encoder if no other TCM encoder 9 = [G,L], for all G G 
Gk.m,!^ and L G Cm, has a superior DS compared to 9. 

An ODS-TCM encoder in Definition |5] is the asymptotically 
optimal TCM encoder in terms of BER and FER for a 
given block length N^. Unlike the more classical definition of 
optimal encoders, ODS-TCM encoders are defined as encoders 
that are optimal in terms of both Aip and B^2. This implies 
that in principle, for some combinations of k,m,v, it is 
possible that no ODS-TCM encoder exists. As we will see in 
Sec.[V] this is not an uncommon situation. Moreover, by using 
this somehow nonstandard definition we avoid listing encoders 
that have optimal BER performance but possibly rather poor 
FER performance (or vice versa). 
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V. Numerical Results 

In this section we study well-structured one- and two- 
dimensional constellations, i.e., AfPAM and A/PSK con- 
stellations. An AfPAM constellation is defined by X = 
[xi,X2,. ■ . ,xmY with Xq = -(M + 1 - 2q)A £ R, 
q = 1,...,M, and A2 = 3/(A/2 - 1) so that = 1. 
An AfPSK constellation is defined by X = [xi\ X2] . . . ; xm] 
with Xq = [cos(27r(q - 1)/Af),sin(27r(g - 1)/A/)] £ 
and q = I, . . . ,M. 

In the following sections we show results of exhaustive 
searches over Qk,m.u and TZm, and thus, these results should 
be understood as a complete answer to the problem of jointly 
designing the feedforward encoder and the labeling for TCM 
encoders. The ODS-TCM encoders presented are obtained by 
comparing the first five nonzero elements in the spectrum, 
which we numerically calculate using a generalization of 
the algorithm presented in ||29l Sec. 12.4.3]0 On the other 
hand, the bounds used to compare with simulation results 
were calculated using 20 terms. The tabulated results are 
ordered first in terms of the output of the MFLSA, then in 
lexicographic order for the memories vi, . . . ^Vk, and then in 
lexicographic order for the encoder matrices G. This ordering 
becomes relevant when there are multiple TCM encoder with 
identical (and optimal) five-term DS. 

A. ODS-TCM Encoders for MPAM 

AfPAM constellations are symmetric around zero. Because 
of this, two TCM encoders based on an A/PAM constellation, 
the first one using the labeling L — [ci; C2; . . . ; cm-i; cm] 
and the second one using a "reverse" labeling L' — 
[cm\ cm-i \ • • ■ ; C2; Ci], are equivalent for any A/. This result 
implies that the number of binary labelings that give nonequiv- 
alent TCM encoders is A/pj/2. Specifically, for m = 2 and 
m = 3 (4PAM and SPAM), only 2 and 120 labehngs need 
to be evaluated, respectively, instead of 24 and 40320 in an 
exhaustive search, see Table U 

To generate only the nonequivalent labelings for 

AfPAM, the MFLSA in Algorithm 1 can be modified as 
follows. Replace Af on lines 5 and 6 with e{index), where 
the integer function e(g) is defined as Af/2 if g = and Af 
otherwise. This has the effect of only generating labelings in 
which the all-zero label is among the first Af/2 positions (i.e., 
the first 4 columns of Table HI] for SPAM). 

1} 1/2 and APAM: The results obtained for f? = 1/2 
and 4PAM and different values of v are shown in Table Hill The 
table reports the DS as well as the labeling and convolutional 
encoder for the ODS-TCM encoders (shown as [ J^^). For 
V = b, however, no ODS-TCM encoder was found, i.e., there 
is no TCM encoder that is optimal in terms of both A^p and 
f3rf2. Instead, we list the TCM encoder with best A^2 among 
those with optimal B^p (shown as [-J^), or vice versa (shown 
as [Y^)- In this table we also include Ungerboeck's encoder^ 
which we denote by [ J^. When Ungerboeck's labeling (NBC) 

'"Note that if more than five elements are considered different ODS-TCM 
encoders might be found. 

"Ungerboeck did not report results for v = 1, and thus, we do not include 
them in the Tables, i.e., we only show the ODS-TCM encoder for v = 1. 




Es/No [dB] 

Fig. 3. BER/FER bounds in (22) and )24t and simulations for Ungerboeck's 
encoders and the ODS-TCM encoders in Table |Ill] for A^s = 1000, 4PAM, 
R = 1/2 (1 [bit/symbol]), and = 4,6. 



or Ungerboeck's convolutional encoder coincide with [ J^^ or 
[■]^, we use the notation [-JUAB qj. [.]UB^ respectively. The 
results in Table [III] show that no gains in terms of MED are 
obtained and that the NBC is indeed the optimal labeling 
for all memories. The key difference between Ungerboeck's 
design and the ODS-TCM encoders is the better multiplicities 
obtained. To compare the gains obtained by the ODS-TCM 
encoders over Ungerboeck's encoders, we show in Fig. [3] their 
BER/FER for = 4, 6. This figure clearly shows the gains 
obtained by using the ODS-TCM encoders which are visible 
not only at high SNR, but also for low SNR values (see, 
e.g., the FER markers for i' = 6). 

2) R ^ 2/3 and SPAM: The results for f? = 2/3 and 
SPAM are shown in Table |IV] For = 1, 2, 3, 4, 6 the reported 
encoders are in the form [ J^^, while for ly — 5 no ODS-TCM 
was found, and we use the same notation as for 4PAM. Unlike 
for f? = 1/2, the parity-check matrix reported by Ungerboeck 
for R = 2/3 specifies the code but not the encoder. To have a 
fair comparison between Ungerboeck's codes with the ODS- 
TCM encoders, we first listed all the convolutional encoders 
that give Ungerboeck's parity-check matrix and then pick the 
one with optimal B^2 (all of them have the same A^2). These 
are the encoders reported in Table |IV] as [-J^. Even though 
Ungerboeck's encoders in Table |IV] are the best encoders for 
that particular parity-check matrix, they coincide with the [ J^ 
encoders only for one out of six cases (i^ = 5). For all the other 
cases, the ODS-TCM encoders result in a better spectrum. 
Also, unlike for 4PAM, Table|IV]shows that the NBC is not the 
optimal labeling. For example, for = 4, the optimal labeling 
is = [1, 2, 4, 0, 6, 5, 3, 7]^^^, which does not follow the SP 
principle (cf. Definition [TJ. In Fig. H) we show the BER/FER 
results obtained by the ODS-TCM encoders for R = 2/3, 
SPAM, and h' — 4,6. This figure shows the tightness of the 
bounds and again gains over Ungerboeck's encoders. 
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TABLE III 

Distance SPECTRUM of ODS-TCM encoders ([ J^^) and Ungerboeck's encoders ([ J^) for A; = 1 [bit/symbol] and 4PAM (m = 2). The 

NOTATION [■]* AND f-]^ IS USED WHEN NO ODS-TCM ENCODER WAS FOUND. 



V 




G 


Distance Spectrum {d^, A^i , -B^2 } 


1 


[0,1,2,3]*'* 


[3,1]*" 


{4.00,0.50,0.50}, {4.80,0.50,1.00}, {5.60,0.50,1.50}, {6.40,0.50,2.00}, {7.20,0.50,2.50} 


2 


[0.1,2,3]™" 


[5,2]" 
[7,2]*" 


{7.20,1.00,1.00}, {8.00,1.25,2.50}, {8.80,1.75,5.25}, {9.60,2.56,10.25}, {10.40,3.81,19.06} 
{7.20,0.50,0.50}, {8.00,1.25,2.50}, {8.80,1.63,4.88}, {9.60,2.56,10.25}, {10.40,3.78,18.91} 


3 


[0,1,2,3]™" 


[13,4]™" 


{8.00,0.25,0.50}, {8.80,1.00,3.00}, {9.60,1.56,6.25}, {10.40,2.75,9.75}, {11.20,3.14,16.84} 


4 


[0,1,2,3]™" 


[23,4]" 
[23,10]*" 


{8.80,0.63,1.88}, {9.60,0.50,2.00}, {10.40,2.00,6.00}, {11.20,2.02,10.09}, {12.00,2.03,13.22} 
{8.80,0.13,0.38}, {9.60,0.50,2.00}, {10.40,1.88,5.38}, {11.20,2.39,10.34}, {12.00,3.72,21.03} 


5 


[0,1,2,3]™" 


[45,10]"" 
[55,4]* 


{10.40,1.13,1.63}, {11.20,1.52,5.09}, {12.00,2.59.12.16}, {12.80,3.58,22.13}, {13.60,5.29,38.60} 
{10.40,0.75,1.75}, {11.20,2.13,8.75}, {12.00,2.14,10.48}, {12.80,4.47,24.75}, {13.60,5.45,37.01} 


6 


[0,1,2,3]™" 


[103,24]" 
[107,32]*" 


{11.20,2.34,5.91}, {12.80,2.82,22.01}, {14.40,7.60,57.35}, {16.00,31.39,268.35}, {17.60,74.37,779.76} 
{11.20,0.13,0.50}, {12.00, 1.44,5.81}, {12.80,1.41,5.77}, {13.60,1.73,12.58}, {14.40,4.58,31.53} 


7 


[0,1,2,3]™" 


[235,126]" 
[313,126]*" 


{12.80,2.19,8.19}, {14.40,3.05,17.66}, {16.00,10.09,89.43}, {17.60,25.03,231.04}, {19.20,90.45,920.63} 
{12.80,1.46,8.02}, {14.40,4.77,34.60}, {16.00,15.42,130.51}, {17.60,35.60,375.08}, {19.20,103.30,1213.89} 


8 


[0,1,2,3]™" 


[515,362]" 
[677,362]*" 


{13.60,0.53,4.66}, {14.40,1.89,10.79}, {15.20,1.66,14.10}, {16.00,3.81,30.45}, {16.80,6.03,49.34} 
{13.60,0.36,2.05}, {14.40,1.06,6.41}, {15.20,1.47,11.09}, {16.00,3.44,23.69}, {16.80,5.25,41.32} 



TABLE IV 

Distance SPECTRUM of ODS-TCM encoders (f ]*^) and Ungerboeck's encoders ([■]'^) for k = 2 [bit/symbol] and SPAM (m = 3). The 

NOTATION l']^ AND [-J^ IS USED WHEN NO ODS-TCM ENCODER WAS FOUND. 



V 




G 


Distance Spectrum {d^ , A^2 , } 


1 


[1,2,4,0,6,5,3,7]*" 


[1,1,1: 1,3,0]*" 


{0.95,1.13,0.84}, {1.14,1.13,1.69}, {1.33,1.13,2.53}, {1.52,1.13,3.38}, {1.71,1.13,4.22} 


2 


[0,1,2,3,4,5,6,7]"*" 


[1,0,0; 0,5,2]" 
[1,0,0; 0,7,2]*" 


{1.71,2.25,1.88}, {1.90,3.52,5.11}, {2.10,6.05,12.35}, {2.29,10.56,27.64}, {2.48,18.47,58.91} 
{1.71,1.69,1.69}, {1.90,3.52,5.11}, {2.10,6.01,12.34}, {2.29,10.56,27.64}, {2.48,18.46,58.91} 


3 


[0,1,2,3,4,5,6,7]" 
[1,2,4,0,6,5,3,7]*" 


[1,0,0; 0,13,4]" 
[1,1,1; 2,15,0]*" 


{1.90,1.27,2.11}, {2.10,3.38,6.75}, {2.29,5.49,14.14}, {2.48,12.45,32.48}, {2.67,18.59,64.81} 
{1.90,1.27,1.90}, {2.10,3.38,8.44}, {2.29,5.49,17.25}. {2.48,12.45,38.50}, {2.67,18.59,74.81} 


4 


[0,1,2,3,4,5,6,7]" 
[1,2,4,0,6,5,3,7]*" 


[1,0,0; 0,23,4]" 
[1,1,1; 2,31,0]*" 


{2.10,2.64,5.59}, {2.29,2.53,6.75}, {2.48,6.75,13.50}, {2.67,12.11,40.55}, {2.86,15.99,66.51} 
{2.10,0.95,1.90}, {2.29,2.53,7.59}, {2.48,7.91,21.78}, {2.67,13.21,45.70}, {2.86,19.77,88.01} 


5 


[0,1,2,3,4,5,6,7]"*" 


[1,0,0; 0,45,10]"" 
[1,0,0; 0,55,4]* 


{2.48,4.32,6.54}, {2.67,7.99,19.45}, {2.86,14.26,46.29}, {3.05,27.05,102.83}, {3.24,44.27,201.33} 
{2.48,3.80,6.96}, {2.67,8.74,21.63}, {2.86,13.53,45.10}, {3.05,29.51,106.50}, {3.24,44.49,198.08} 


6 


[0,1,2,3,4,5,6,7]"*" 


[1,0,0; 0,103,24]" 
[1,0,0; 0,107,32]*" 


{2.67,10.74,22.97}, {3.05,19.91,86.93}, {3.43,72.68,343.40}, {3.81,353.99,1927.40}, {4.19,1137.86,7442.94} 
{2.67,1.42,4.27}, {2.86,8.46,24.43}, {3.05,12.94,40.47}, {3.24,15.68,74.20}, {3.43,40.61,182.47} 



B. ODS-TCM Encoders for MPSK 

A TCM encoder based on an AfPSK constellation is not 
affected by a circular rotation of its labeling, i.e., without 
loss of generality it can be assumed that the all zero label 
is assigned to the constellation point Xi = [1,0]. The conse- 
quence of this is that for A/PSK constellations, the number 
of reduced column echelon matrices that give nonequivalent 
TCM encoders is further reduced by a factor of M. In view 
of the results in Table U for 4PSK, there is only one labeling 
that needs to be tested, e.g., the NBC. For m > 3, the 
nonequivalent labelings can be obtained from the MFLSA by 
setting index <— 3 in line 4, which gives the FLSA of ll28l . 
For example, for M = 8, the output corresponds to the first 
column of Table [III which gives 30 labelings. 

1} R = 1/2 and APSK: In this case there is only one 
labeling to be tested (the NBC), and thus, only a search over 
the encoders needs to be performed. Moreover, without loss 
of generality, we can use the BRGC instead (because it is 
in the same Hadamard class as the NBC) and search over 
encoders for this labeling. Since 4PSK with the BRGC can 
be considered as two independent 2PAM constellations (one 
in each dimension), the design of TCM encoders in this case 



boils down to selecting convolutional encoders with optimal 
spectrum (in the sense of Definition |5]l. 

We have performed an exhaustive search for convolutional 
encoders with optimal spectrum up to = 12 and found that 
our results coincide with those reported in [38. Table I]. For 
V ~ 1,2, 3, 4, 5, 6, 11, 12 the optimal convolutional encoders 
([•]^^) are in fact the encoders from ITS', Table I] (which were 
initially optimized only in terms of B^2). For v — 7,8,9, 10 
we found that no optimal encoder exists, i.e., the convolutional 
encoders optimal in terms of ^1^2 are not optimal in terms of 
_B^2 and vice versal3 These encoders are in fact shown in [38l 
Table iQ which extends the results in lfT4l - lfT6l because it 
considers both Aj2 and Bj2 as optimization criteria. 

Based on the discussion above, we conclude that an ODS- 
TCM encoders can be constructed by concatenating the en- 
coders in ll38l Table I] with a 4PSK constellation labeled 
by the BRGC. Alternatively, ODS-TCM encoders can be 
obtained by using a 4PSK constellation labeled by the NBC 

'^Convolutional encoders with optimal A^2 and memories up to = 26 
have been recently published in t39l Table 7.1]. 

"Although the search in |38| was performed only considering events at 
minimum Hamming distance and not over the whole spectrum. 
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Fig. 4. BER/FER bounds in )22t and )24t and simulations for Ungerboeck's Fig. 5. DS for encoders with u = 4 for R = 2/3 and 8PSK from Table Fvl 
encoders and the ODS-TCM encoders in Table |IV] for A^s = 1000, SPAM, 
= 2/3 (2 [bit/symbol]), and = 4, 6. 



and using the encoders in ll38l Table I] after applying the 
transformation T^^ = [1,1; 0,1]. For example, for v = 8, 
we found G[5i5.677] and G[435_657] to be the optimal encoders 
in terms of and Bj^2, respectively, and thus, the two pairs 
of equivalent ODS-TCM encoders are 8 = [G[5i5,677] j -B2] 
and 6 = [G[5i5_677]T'~\ -^2], and 9 = [G[435_657] , -B2] and 

= [G[435,657]T'^ ,-^2]- 

2) R = 2/3 and 8PSK: The results obtained for i? = 2/3 
and 8PSK are shown in Table |V] Somehow disappointingly, 
this table shows that the NBC is indeed the optimal labeling 
in all the cases, and thus, the selection of the labeling for 
this particular configuration does not provide any gains over 
Ungerboeck's TCM schemes. The better spectrum obtained by 
the ODS-TCM encoders in this case then comes only from the 
selection of the convolutional encoder. 

In Fig.|5] we show the DS for the encoders in Table [Vl with 
v = A. It is clear from the figure that an encoder optimal in 
terms of Aii2 can be suboptimal in terms of Bi^-i, and vice 
versa. In addition, the figure shows how the set of SEDs V 
is in general different for different encoders. It also shows 
how Ungerboeck's encoder is optimal in terms of A112 for the 
term at MED, but in general suboptimal if the whole DS is 
considered. 

We note that depending on the ODS-TCM encoders in 
Table IV] have inferior, equivalent, or superior Bci2 spectrum to 
those listed in |22, Table 3.2], [19, Table 6.10]. The reason 
for this is that the codes tabulated in ^ Table 3.2], fl9\ 
Table 6.10] are found by searching over parity check matrices 
and then converted to feedback encoders (in observer canoni- 
cal form |fT9l Fig. 2.2]). On the other hand, we search over a 
different set of encoders, namely, over all the noncatastrophic 
feedforward encoders. 

All labelings we found for the ODS-TCM encoders (i.e., the 
highlighted labelings in Table |ll] and the optimal ones in 



Tables [TVl and IVT) have optimal EP. This makes us conjecture 
that good TCM encoders can be found by using the EP of |27l 
on top of the proposed classification. This approach would 
indeed reduce the search space (for example, for SPAM and 
8PSK constellations, only eight and two labelings, respec- 
tively, would need to be tested). However, it would not allow 
us to claim optimality in the sense of Definition |5] 

VI. Conclusions 

In this paper we analyzed the problem of jointly designing 
the feedforward convolutional encoder and the labeling of a 
TCM encoder. It was shown that the number of labelings that 
need to be checked can be reduced if they are grouped into 
modified Hadamard classes. This classification allowed us to 
prove that it is always possible to design a TCM encoder based 
on the BRGC with identical performance to the one proposed 
by Ungerboeck in 1982. The numerical results show that in 
most cases, the NBC is the optimal binary labeling for TCM 
encoders and that gains up to 0.3 dB over the previously best 
known TCM schemes can indeed be obtained. 

The classification of labelings presented this paper does 
not make any assumption on the channel nor on the receiver 
Because of this, the presented design methodology can be used 
to design optimal TCM encoders for other channels as well 
as for suboptimal (BICM) decoders. 

The algorithm introduced in this paper to find all the label- 
ings that need to be tested in an exhaustive search becomes 
impractical for constellations with more than 16 points. In this 
case, a suboptimal solution based on an algorithm (inspired by 
the linearity increasing swap algorithm of 1281 Sec. IX]) that 
generates a subset of (good) labelings could be devised. This 
approach could also be combined with the concept of labelings 
with optimal EP |27|. The design of such an algorithm is left 
for further investigation. 
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TABLE V 

Distance SPECTRUM of ODS-TCM encoders ([ J^^) and Ungerboeck's encoders ([ J^) for k = 2 [bit/symbol] and 8PSK (m = 3). The 

NOTATION l]^ AND l]^ IS USED WHEN NO ODS-TCM ENCODER WAS FOUND. 







G 


Distance Spectrum {d^ , A^2 , B^2 } 


1 


[0,1,2,3,4,5,6,7]*'' 


[1,0,0; 0,1,2]*" 


{2.59,2.00,1.50}, {3.17,2.00,3.00}, {3.76,2.00,4.50}, {4.00,1.00,0.50}, {4.34,2.00,6.00} 


2 


[0,1,2,3,4,5,6,7]™"* 


[1,0,0; 0,5,2]™ 


{4.00,1.00,0.50}, {4.59,4.00,4.00}, {5.17,8.00,14.00}, {5.76,16.00,38.00}, {6.34,32.00,96.00} 


3 


[0,1,2,3,4,5,6,7]™ 


[1,2,0; 4,1,2]" 
[1,2,0; 4,5,2]*"* 


{4.59,2.00,2.50}, {5.17,4.00,8.50}, {5.76,8.00,25.00}, {6.00,1.00,0.50}, {6.34,16.00,66.00} 
{4.59,2.00,2.00}, {5.17,4.00,8.50}, {5.76,8.00,25.00}, {6.00,1.00,0.50}, {6.34,16.00,66.00} 


4 


[0,1,2,3,4,5,6,7]™ 


[2,7,0; 7,3,2]"^ 
[2,7,0; 7,1,2]* 
[1,4,2; 6,1,0]® 


{5.17,2.25,5.50}, {5.76,4.63,14.13}, {6.00,1.00,0.50}, {6.34,6.06,26.50}, {6.59,4.00,5.50} 
{5.17,2.25,5.00}, {5.76,3.88,11.56}, {6.00,1.00,0.50}, {6.34,9.56,38.81}, {6.59,4.00,5.50} 
{5.17,2.50,5.00}, {5.76,3.75,11.25}, {6.34,8.13,32.44}, {6.59,3.50,4.50}, {6.93,16.19,80.94} 


5 


[0,1,2,3,4,5,6,7]™ 


[1,2,0; 30,25,16]" 
[1,2,0; 30,25,10]*"* 


{5.76,4.00,10.50}, {6.00,1.00,0.50}, {6.34,4.00,16.25}, {6.93,4.00,24.13}, {7.17,3.00,7.50} 
{5.76,2.00,5.75}, {6.00,1.00,0.50}, {6.34,3.63,15.56}, {6.59,3.00,5.50}, {6.93,8.06,40.63} 


6 


[0,1,2,3,4,5,6,7]™ 


[4,11,0; 13,4,6]" 
[1,6,0; 27,25,12]* 
[1,6,0; 35,31,6]"* 


{6.34,5.25,22.56}, {7.17,10.00,28.88}, {7.51,14.53,98.50}, {8.00,3.00,3.75}, {8.34,38.56,199.78} 
{6.34,3.25,12.00}, {7.17,7.25,17.88}, {7.51,19.13,119.17}, {8.00,3.00,5.00}, {8.34,36.69,159.69} 
{6.34,3.56,11.50}, {7.17,7.25,16.88}, {7.51,16.58,92.05}, {8.00,3.50,4.75}, {8.34,30.63,150.81} 
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